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Within a two-band model for the recently discovered ferropnictide materials, we calculate the 
thermal conductivity assuming general superconducting states of A\ g ( "s-wave" ) symmetry, consid- 
ering both currently popular isotropic "sign- changing" s states and states with strong anisotropy, 
including those which manifest nodes or deep minima of the order parameter. We consider both 
intra- and interband disorder scattering effects, and show that in situations where a low-temperature 
linear-T exists in the thermal conductivity, it is not always "universal" as in d-wave superconduc- 
tors. We discuss the conditions under which such a term can disappear, as well as how it can be 
induced by a magnetic field. We compare our results to several recent experiments. 

PACS numbers: 



INTRODUCTION 



The symmetry class of the newly discovered ferrop- 
nictide superconductors 2 - is still controversial, due in 
part to differing results on supcrfluid density^ 
angle-resolved photoemission (ARPES) 
clear magnetic resonance (NMR)^^^^^^, Andreev 
and other experimental probes. 
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spectroscopy 

In some cases, these experiments have been interpreted 
as implying the absence of low-energy excitations, i.e. a 
fully developed spectral gap. In others, low-energy ex- 
citations have been observed, and taken as indication of 
the existence of order parameter nodes. It may be that 
these differences depend on the stoichiometry or doping 
of the materials, which affects the pairing interaction, or 
sample quality, or both. 

As in other classes of potentially unconventional super- 
conductors, one's ability to identify the symmetry class of 
a candidate material by observation of low-T power laws 
in temperature reflecting low energy quasiparticle exci- 
tations is limited by how low in T one can measure. At 
intermediate temperatures, variations of thermodynamic 
and transport properties can be affected by details of 
band structure, elastic and inelastic scattering, as well as 
the presence of thermal phonons. Only at the very low- 
est T can one - in principle - extract direct information 
on the order parameter structure. Thermal conductiv- 
ity measurements have played an important role in past 
discussions of unconventional superconductivit y 25 ' 26 , in 
part because they can be extended to T of order tens of 
mK. In addition, such measurements are distinguished 
because they are bulk probes, and because they are 
unusually sensitive to the presence of order parameter 
nodes. If lines of nodes are present, the thermal con- 
ductivity k(T) manifests a low-T linear, in temperature, 
term which is purely electronic in origin and is associated 
with residual quasiparticle states at the Fermi level, in- 
duced by disorder or a magnetic field. If, in addition, the 
order parameter averages to zero over the Fermi surface 



(as in the d-wave case appropriate for the cuprates), this 
linear-T term in zero field is known to be "universal" , 
in the sense that its magnitude is only weakly disorder 
dependent. 

Very recently, several low-temperature measurements 
of thermal transport have been made on the BaFe2As2 
(Ba-122) material doped with C<M and Ni 2 ^, as 

well as on the stoichiometric superconductor LaFcPO 3 -. 
In the case of the Ba-122 samples, cither zero or very 
small lincar-T terms have been reported in zero field, 
leading to the conclusion that there is a fully developed 
spectral gap in these materials 3 ^ in these experimental 
works. This is in contrast to reports of power law tem- 
perature dependence in the superfhud density measured 
on the same materials 3 ^! 5 -^, as well as other strong in- 
dications of low-energy excitations. One way to reconcile 
these experiments is to note that thermal conductivity at 
mK temperatures probes lower energy scales than those 
measured in other experiments to date; thus it is possi- 
ble that a band of low-energy excitations extends to very 
low energies, but not all the way to the Fermi level, ei- 
ther due to an intrinsic highly anisotropic order parame- 
ter with deep minima, or a band of impurity states which 
lies at low but nonzero excitation energies. Such impurity 
states can be produced, e.g., in isotropic, sign-changing s- 
wave pair stat o 33 : 34 ' 35 ' 36 allowable in multiband systems 
if special conditions on the ratio of intra- to interband 
scattering are met. However, a further strong constraint 
from the thermal conductivity measurements is that a 
significant linear-T term in the thermal conductivity is 
observed with the application of a small magnetic field 
of order one Tesla and hence much below the upper crit- 
ical field, H C 2- This residual term grows with increas- 
ing field. This would be consistent with the existence 
of quasiparticle states at low, but finite energy. Refs. 
12811371 on K- and Co-doped Ba-122 samples claimed that 
this enhancement is significantly larger than that to be 
expected in the case of a conventional s-wave supercon- 
ductor. In contrast, Ref. l29lcame to opposite conclusions 
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on the measurements of a Ni-doped sample, and Ref. 
reported a small but significant linear-T term in zero field 
in Co-doped Ba-122. 

A recent measurement of k(T) on the fcrrophos- 
phide superconductor LaFcPO finds a very large linear-T 
term2i. If this is interpreted as indicative of order pa- 
rameter nodes, it would be consistent with the linear-T 
dependence in the superfluid density also observed for 
this system^. Note that disorder in a sign-changing s 
state cannot produce such a term in the superfluid den- 
sity. This material is the only material yet discovered 
among the growing ferropnictide family of superconduc- 
tors whose undoped "parent compound" is superconduct- 
ing at zero pressure. It is therefore expected to be signifi- 
cantly cleaner than other superconductors discussed here. 
This may be relevant because it has been proposed that 
disorder in highly anisotropic "s-wave" (Ai ff symmetry) 
states can "lift" shallow nodes in the order parameter, 
leading to a fully developed spectral gap^. The authors 
of Ref. [HI note that, despite a very sharp resistive tran- 
sition, the cannot completely exclude the possibility that 
the lincar-T term is partly extrinsic; however, even in 
that case the dominant dependence of the thermal con- 
ductivity on the magnetic field should come from the 
superconducting phase. 

There is a developing consensus that the gap changes 
sign between the electron and hole Fermi surface sheets. 
From the theoretical standpoint, states with nodes or 
deep minima appear to be quite natural. Several micro- 
scopic theories of the spin fluctuation mediated pairing 
interaction in the ferropnictide materials have attempted 
to calculate the momentum dependence of the order pa- 
rameter associated with the leading superconducting in- 
stability. Using a 5 Fe-orbital model, Kuroki et al^ 
performed an RPA calculation of the interaction to con- 
struct a linearized gap equation, and determined that 
the leading pairing instability had s-wave {Ai g ) symme- 
try, with nodes on the electron- like Fermi surface ("/3 
sheets"). Wang et al£2- studied the same pairing problem 
within a 5-orbital framework using the functional renor- 
malization group approach, also finding that the leading 
pairing instability is in the Ai 5 -wave channel, and that 
the next leading channel had B\ g (d x i_ y i) symmetry. 
For their interaction parameters, they found no nodes 
on the Fermi surface, but nevertheless a significant vari- 
ation of the magnitude of the gap. Graser et al also 
performed a 5-orbital RPA framework^, using the DFT 
bandstructure of Cao et al4£ as a starting point. These 
results indicated that the leading pairing channels were 
indeed of s (Ai g ) and d x 2_ y 2 symmetry, and that one or 
the other could be the leading eigenvalue, depending on 
details of interaction parameters. More recently, several 
authors have investigated the factors including intrasheet 
Coulomb interaction, nesting of electron pockets, and or- 
bital character of pairing which can influence order pa- 
rameter anisotropy within these model o 42 ! 43 ' 44 . Other 
approaches have also obtained A\ g gaps which change 
sign between the hole and electron Fermi surface sheets 



but remain approximately isotropic on each shee t ' . 

In this paper we calculate the expected thermal con- 
ductivity in superconducting states potentially appropri- 
ate to the ferropnictide superconductors. We adopt for 
convenience a phenomenological 2-band model, allowing 
order parameters on two Fermi surface sheets represent- 
ing the hole- and electron- doped sheets found in density 
functional theory calculations for these materials. We 
consider the region outside of the doping range where 
superconductivity may coexist with antiferromagnetism. 
Our model for disorder consists of terms allowing for scat- 
tering within (intraband) and between (interband) Fermi 
surface sheets, of arbitrary strength. This allows us to 
control the width and position of the impurity band in 
both nodal pairing states and those with a fully devel- 
oped spectral gap, which we examine with a view towards 
determining the size and universality of the linear-T term 
in k at the lowest temperatures. After examining the 
zero-field situation, we discuss the effect of an applied 
field. To this end we adopt the method of Brandt-Pesch- 
Tewordt to obtain predictions for the widest possible field 
range. We illustrate the various possibilities of supercon- 
ducting state and disorder types which allow the results 
observed thus far in experiments. 

II. MODEL 

We begin by assuming a metallic system with two 
bands 1 and 2, characterized by densities of states Ni 
and N2 at the Fermi level, and a pair interaction which 
is a sum of separable terms, 

V(k,k') = y 1 $ 1 (k)$ 1 (k') + V r 2 $ 2 (k)$ 2 (k') 

+ V 12 [$!(k)$ 2 (k') + $ 2 (k)$i(k , )] ) (1) 

where $i is function of A\ g symmetry depending on mo- 
mentum restricted to band i = 1,2. 

For disorder we will assume an orbital-independent 
matrix element which scatters quasiparticles either 
within a given band with amplitude Ua, i = 1,2, or be- 
tween bands with amplitude U12. As discussed in the 
appendix, we sum all single-site scattering processes of 
arbitrary strength to obtain a disorder-averaged Nambu 
self energy E = ni mp T, where ni mp is the concentration 
of impurities. For simplicity, we assume Un = £/ 22 = 
Ud, with equal densities of states Ni = N through- 
out the paper. In our preliminary considerations we re- 
strict ourselves to purely intraband scattering, Uyi = 0. 
The disorder is characterized by two intraband scatter- 
ing parameters on each sheet: Ti = rii mp / (jr N) and 
d = l/(irNiUu); For our simple initial case with 2 sym- 
metric bands we set T; = T and Ci = c, i = 1,2. The 
initial neglect of interband scattering may be understood 
in zeroth order by noting that a screened Coulomb po- 
tential with screening length of order a unit cell size 
will gencrically have larger small-q compared to large-q 
scattering. The real situation is somewhat more com- 
plex since the same orbitals contribute to both electron 
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and hole Fermi surface sheets^, and therefore a sub- 
stitutional impurity, such as Co, may be expected to 
produce a significant interband scattering component as 
well. Hence we relax this requirement and below also an- 
alyze the regime U\2 — Ud, a nd, in particular, the case 
U12 = Ud, as discussed in Refl35l. Weak interband scat- 
tering U12 -C fii, U22 does not qualitatively change the 
results obtained in the limit U12 = 0. 

The full matrix Green's function in the presence of 
scattering in the superconducting state is given by a diag- 
onal matrix in band space, as discussed in the appendix, 



G(k,w) 



grp + e k T 3 + A k ri 

~,2 ?2 A 2 : 



(2) 



where k = 6 Si is restricted to Fermi surface sheet Si 
with i = 1,2, and the renormalized quantities Co = lo — Y,q, 
£k = £k + S3, Ak = Ak + Si also depend on the band 
indices through k. The S Q are the components of the self- 
energy proportional to the Pauli matrices r Q in particle- 
hole (Nambu) space. 

Below we focus on two quantities. It is useful to start 
with the analysis of the total density of quasiparticle 
states (DOS) 



N(co) 



2 ^TrIm^G(k,w) 

k 



(3) 



i ki 



where the second equality indicates the explicit integra- 
tion over distinct Fermi surface sheets with momenta k^ . 
We will be comparing results for the DOS with thermal 
conductivity k calculated using the standard approach^, 



V NjV 2 Fi 



LO n ^ 

dw^sech (— ) 



1 



I A* 



(4) 



Here < .. >^ is average over the Fermi surface sheet, 
denotes the sum over the bands, and VFi is the Fermi 
velocity on sheet i. We assume cylindrical Fermi surfaces, 
so that v Fi is isotropic. 



III. ISOTROPIC Ai 9 STATES 



We first discuss the thermal conductivity in isotropic 
s-wave or Ai g states. If the sign of the order parameter 
A is the same on both sheets, and no magnetic disorder 
is present, the low-T thermal conductivity will be similar 
to classic calculations for conventional superconductors, 
and yield an exponential low-T dependence for the elec- 
tronic part. In the case of a sign-changing s ("s±") state 
proposed by Mazin et al4^, the situation is more inter- 
esting. Here one assumes an isotropic A, on each sheet 
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FIG. 1: Density of states (a) and thermal conductivity (b) for 
an isotropic s± state with Ai = — A2, shown for Ud = Un = 
U22 (intraband Scattering ) and scattering rate parameters 
c = 0.07 and T = 0.3T c o in cases (i) weak intraband scattering 
only, Ui2/Ud=0 (solid line); (ii) pairbreaking scattering with 
midgap impurity band, Un/Ud = 0.98 (dashed line); (iii) 
pairbreaking scattering with impurity band overlapping Fermi 
level, Ui 2 /U d = 1.0 (dotted line). 



i, but assumes that sgnAi = — sgnA2. In terms of Eq. 
[2 we choose the functions $;(k) = 1 for k <E Si and zero 
otherwise, and fix the sign of V12 to be opposite to that 
of V\ , V2 so as to induce a sign-changing order parameter 
between the two sheets. In the clean case, we continue 
to expect an exponential or full gapped dependence to 
the thermal conductivity. On the other hand, in such 
systems ordinary disorder is pairbreaking if it includes 
a strong interband componen t 33 ' 34 i 35 ' 36 . For some situ- 
ations a low-energy impurity band may indeed give rise 
to a linear term in the thermal conductivity. 

For simplicity, we assume that A2 = — Ai = A, and 
equal densities of states on the two bands. In Fig. [U 
we now illustrate the correspondence between the for- 
mation of the impurity band in the fully gapped state, 
and the creation of the linear term. In the absence 
of interband scattering, there is no pairbreaking in the 
sense of Anderson's theorem, and the spectral gap in the 
DOS is identical to the unrenormalized order parameter 
A, corresponding to an activated thermal conductivity 
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FIG. 2: Density of states (insert) and thermal conductivity 
for an isotropic s± state with scattering parameters as Fig. Q] 
but U 12 /U d =1.0 and T = 0.2,0.25, and 0.3T c0 . 



~ exp(— A/T). As interband scattering is increased, 
states are pulled down from the continuum into the gap, 
creating eventually a band of midgap states in the DOS, 
as shown. If there is still a narrow energy range which is 
gapped near the Fermi level, this smaller gap determines 
the slower but still exponential decrease of k at the lowest 
temperatures. As soon as the impurity band of midgap 
states overlaps the Fermi level, a linear term in k ap- 
pear a 34 ' 55 . Note that a significant interband component 
of scattering essentially equal to the intraband compo- 
nent is absolutely necessary for this to occur, which re- 
quires special conditions as described above. 

We note further that if one varies the concentration of 
impurities in the situation with U12 = Ud, the change in 
the residual density of states N(0) is reflected directly in 
the slope of the thermal conductivity, as shown in Fig. [5] 
This is not surprising but is in dramatic contrast to the 
universal (disorder-independent) behavior observed in d- 
wave superconductors. This raises the question of the 
degree of universality of transport coefficients in pairing 
states which have s-wave (A lg ) character but also display 
nodes. 



IV. 



ANISOTROPIC Ai s STATES 



We now examine states within the same A lg symmetry 
class, but where gap minima arc cither very deep, with no 
sign change, or with actual sign change (nodes). To make 
contact with microscopic theory (see e.g. Ref. Ell ), we 
assume that one of the sheets (in microscopic theory, the 
so-called "a" sheet around the T point) has an isotropic 
order parameter while the other (the "/?" sheet around 
the M point) has a highly anisotropic one. In this case, 
the order parameters in the two bands are given by 



Ai 



A 2 = A iso + A am cos2</>, 



(5) 
(6) 



where <fi is the angle around the electron Fermi surface 
sheet 2. It useful to define the gap ratio in the electron- 
like band, r = A ara j/Aj SO , so that r > 1 (A Qn j 

> A lso ) 

gives a state with nodes in that band, while r < 1 
(A ara < A iso ) has none. 

To obtain this order parameter from Eq. (JTJ) , we choose 



*i(k) = 
$ 2 (k) = 



1 keSi 

otherwise 

1 + r v cos(2<^) ke5 2 

otherwise 



(7) 

(8) 



where as before, S± and S2 represent the hole and electron 
Fermi surfaces, respectively. In the clean limit, the gap 
equation is 



A,(0) = 2irTj2M0) 



:(9) 



where w„ are fcrmionic Matsubara frequencies. Below 
we adjust the value of rv to study a nodal system with 
r = 1.3 and an anisotropic state with no nodes with 
r = 0.9. 

In the presence of disorder, we evaluate the impurity 
average self energies E i a for both the bands, where again 
i is the band index and a is the Nambu index. This 
calculation is detailed in the Appendix, with the results 
presented in Eqs. (|A8[ )- (|All[) . Since the structure of the 
order parameter in the clean limit already supports low- 
energy excitations, we first ignore the interband scat- 
tering in the first analysis, and focus on the effects of 
intraband scattering alone. We define the conventional 
renormalized quantities 



A, 



w - S ij0 i 
A, + E u , 



1,2 



(10) 
(11) 



where in each case the first subscript is a band index, 
while the second one is a Nambu index. Since the self 
energy is k-independent in this approximation, we can 
associate the self-energy with the renormalization of the 
isotropic component, Ai SO — Ai SO + £24, but this is sim- 
ply a matter of convenience. The total thermal conduc- 
tivity comes from sum over both bands, but at very low 
temperatures the contribution from the first band is very 
small due to the fully developed gap assumed. 
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states with nodes 



We first discuss the situation where A 2 (k) has nodes 
but a non-zero average over the Fermi surface, and for 
concreteness take r = 1.3. The behavior of the low-T 
thermal conductivity with increasing disorder for a case 
with individual scatterers near the strong potential limit 
is now shown in Fig. [3] 
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FIG. 3: Density of states N(w)/NTotai (top) and normal- 
ized thermal conductivity (k(T)/T)/(k„/T c ) vs. T/T c for 
two band anisotropic model with isotropic order parameter 
on sheet 1 and anisotropic order parameter Ai = —1.1 T c o 
on sheet 1 and Ai SO = 1.3 T c o , r = 1.3 on sheet 2 (Eq. JSJ). 
Results shown for various values of intraband scattering rate 
T/Tco and c = 0.07 , and no intraband scattering, U12 = 0. 



The evolution of k(T) is very different from that for 
a pure ri-wave superconductor. This is clearly seen from 
the evolution of the T = limit of the thermal conduc- 
tivity. In the pure case the linear term is nearly invisible. 
As intraband disorder is increased, the linear term signif- 
icantly increases in magnitude, goes through a maximum 
and eventually disappears, leading to an exponential tem- 
perature dependence. To some extent this behavior can 
be understood by examining the corresponding density 
of states, as shown in the upper panel; as disorder in- 
creases, the nodal quasiparticle states are broadened and 
a residual density of states appears, but as disorder is in- 
creased further the nodes are lifted and a fully developed 
spectral gap appears, as discussed in Ref. l38l 

To clarify why there is no "universal indepen- 
dence" of weak disorder expected, e.g. for d wave 
superconductor o 49 ' 50 ! 51 , we plot in Fig. [4] the value of 
the asymptotic low-T limiting value of k/T as a function 
of disorder; there is, for this case, no range of disorder 
where the behavior can in any sense be called univer- 
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FIG. 4: Magnitude of linear-T term in thermal conductivity 
in limit T — > for case with nodes r — 1.3, plotted as a 
function of intraband scattering rate T/T c o for c = 0.07. Solid 
line: exact numerical result. Dashed line: analytical estimate 
from text, Eq. I)14p . Here koo is the thermal conductivity 
for a d— wave superconductor, with order parameter A(<^>) = 
A a „i cos 2(f>. Aani is same as anisotropic component in clean 
system for the model specified by Eq. ([6]). 



sal. This result is somewhat similar to that in Ref. [52] 
where the effect of an orthorhombic distortion on the in 
plane thermal conductivity of YBCO was studied. Note 
that in their case the s— wave component of the d + s or- 
der parameter breaks A\ g symmetry in the ab-plane. In 
our case A\ g is preserved because there are two (3 sheets 
whose nodal structures are rotated with respect to each 
other by 90 degrees 4 ^. 

To analyze the origin of the breakdown of universal- 
ity in the anisotropic A\ g state, we evaluate the T — » 
limit by replacing the derivative of the Fermi function 
by a delta function and integrating over the energy uj. 
The isotropic band 1 does not contribute to the thermal 
conductivity at very low temperature because it is fully 
gapped. The main contribution comes therefore from the 
nodal states from band 2. In contrast to a ri-wave super- 
conductor, here the anomalous self-energy £2,1 is finite, 
therefore the compensation between the density of states 
and the scattering rate does not occur, and the universal 
behavior breaks down. As a result, the position of the 
nodes on the Fermi surface shifts, and the slope of the 
gap changes. 

If we linearize the gap near the node, A(0) w 
&f^a(0o ~ < t ) )i where r/>o is the location of node, deter- 
mined from cos2^o = -[A iao + £ 2 .i(w = 0)]/A am , the 
renormalized gap slope is 



DA = 2k p 1 A ani sin(2</>o 



(12) 



2kp 



(A iso + S 2 , 1 (a J = 0)) 2 . (13) 



6 



Summing over the nodes, we find 



K 
f 



Nv 2 f tt 2 

3 kFVA 



(14) 



which has precisely the same form as the well-known d- 
wave result (to leading order in v^/vf), except that the 
gap velocity va, which is unrenormalized by disorder in 
the d-wave case, is strongly disorder dependent here due 
to the nonzero off-diagonal impurity self-energy £2.1 in 
Eq. [T31 The increase in the residual thermal conduc- 
tivity is therefore due to the flattening of the gap in the 
near nodal region before the system becomes fully gapped 
at higher impurity scattering rates. The absence of the 
residual linear term in this picture is only consistent with 
a sufficiently high disorder, when the spectral gap is fi- 
nite. 



B. Anisotropic states with deep gap minima 

An alternative scenario for the absence of a linear 
term in k(T) in K-, Ni- and Co-doped 122 fcrropnictidc 
material s 28 ^ 9 ' 37 , is a highly anisotropic state on at least 
one of the Fermi surface sheets with deep minima but no 
true nodes. Here we choose A an i < Ai SO (r = 0.9) to sim- 
ulate a situation where the clean state is slightly gapped, 
and varies between a minimum value of Ai SO — A an i and 
A.; so + A an i . Again we begin by including only intraband 
impurity scattering, for which the results are shown in 
Fig. In this case, disorder merely increases the spec- 
tral gap due to averaging, as discussed in Ref. l38l . leading 
to an increasingly rapid exponential decay 

If interband scattering is included, low-energy states 
appear, similarly to the s± case considered above, but 
again we find that the values of interband scattering 
strength U12 near the intraband value Ud are necessary 
to create such states sufficiently near the Fermi level to 
create a linear term in k. For simplicity, therefore, we 
take both inter and intra band scattering potentials to be 
equal, i.e. the scattering is isotropic in momentum space. 
For intermediate to strong potentials, states are then cre- 
ated near the Fermi level. It is important to remember 
that such scattering rapidly suppresses T c , as illustrated 
in Fig. [5] The next figure exhibits the thermal con- 
ductivity for this system. Because the pure system has a 
small spectral gap, even the smallest disorder in this limit 
gives rise to an impurity band close to the Fermi level, 
creating unpaired quasiparticles. For significant impu- 
rity concentrations, a strong linear term appears which 
also of course violates universality, as shown explicitly in 
Fig. [Jj The variation of this linear term with disorder 
for isotropic scattering is shown in Fig. [H] and compared 
to an approximation where we made a series expansion 
around ir/2 for A^. We find 




FIG. 5: (Top) Normalized density of states N(w)/No vs. 
a)/T c o for two band anisotropic model with isotropic gap on 
sheet 1 and r — 0.9 (deep gap minima) on sheet 2 (Eq. JB). 
Results are shown for various values of intraband scattering 
rate T/T c o and c = 0.07. Bottom: normalized thermal con- 
ductivity (k(T)/T)/(k„/T c ) vs. T/T c . 
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1 
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6 (r 1 2 + (A + s 1 , 1 ) 2 ) 3 / 2 ' 



(15) 



where Ti are the normal state scattering rates defined 
above, and all the self-energies are evaluated at w = 0. 
Note that there now appears a contribution from the 
isotropic band 1 because Eq. [T5] assumes the isotropic 
scattering condition Un = U12, which leads to strong 
pairbreaking and quasiparticle states near the Fermi 
level. Consequently, the absence of the residual linear 
term in k is also consistent with the deep minima pro- 
vided the interband scattering is not too strong. We now 
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FIG. 6: Comparison of the effects of intraband and isotropic 
scattering on T c . Solid curve is pure intraband scattering 
Tc/T c o vs. r/Tco for U12 = 0; dashed curve is same, but for 
U12 = Un = Uii- 



proceed to investigate the field dependence in the two 
cases. 



V. FIELD DEPENDENCE 

Thermal conductivity depends on the applied mag- 
netic field since the density of unpaired electrons depends 
on the field magnitude. These electrons carry entropy 
and hence enhance the heat current. They also scatter 
phonons and therefore reduce the lattice contribution to 
the thermal transport, so that the two effects compete. 
On general grounds, n(T,H) that increases at low tem- 
peratures with applied field can be assumed to contain 
a substantial electronic component^. In some systems, 
such as heavy fermion metals, the electron contribution 
to the thermal conductivity is dominant, allowing a direct 
probe of the heat transport in the superconducting state 
throughout the T-H planed. In other materials, where 
the phonon contribution is substantial, the quantity that 
lends itself most easily to analysis is the residual linear 
term in the thermal conductivity, limT_>o K /T, which is 
purely electronic since the phonon contribution vanishes 
in that limil^iS^. 

Therefore, for the purposes of comparison with exper- 
iment, we focus on the field dependence of the electronic 
thermal conductivity at low temperature. In nodal su- 
perconductors, where the transport is dominated by bulk 
quasiparticlcs with momenta nearly along the nodal di- 
rections, two methods have been employed to describe 
this dependence. The semiclassical approach is based on 
the observation that the energy to break a Cooper pair 
is lowered outside of the vortex core since the unpaired 
electrons do not participate in the supcrcurrent flowing 
around the vortex. Hence effect of the field can be de- 
scribed by the Doppler shift of the quasiparticlc energy, 
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FIG. 7: Density of states N(w)/No (top) and normalized 
thermal conductivity k(T)T c /(k„T) vs. T/T c for two band 
anisotropic model with isotropic order parameter on sheet 1 
and r = 0.9 (deep gap minima) on sheet 2 (Eq. ©. Re- 
sults are shown for equal intraband and interband scattering 
Uia = U d . 




FIG. 8: k/T in T — > limit, for the state with deep gap 
minima. «oo is thermal conductivity for pure d-wave state. 
Band 1 has isotropic order parameter A = —1.7 T c o and band 
2 has r — 0.9 and Ai SO = 1.5 T c q on it. 
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v s (r) ■ hp, where v s (r) is the supcrvclocity field deter- 
mined by the vortex structur e) 77 ' 78 ! 79 . This energy shift 
is local, and therefore the method is very well suited for 
describing the thermodynamic quantities, but requires 
additional assumptions to account for correlation func- 
tions and transport propertie d 80 ' 81 i 82 i 83 , It is applicable 
at low energies and therefore restricted to low tempera- 
tures and fields. 

An alternative approach assumes the existence of the 
vortex lattice and describes the behavior starting from 
the moderate to high field regime. The approximation 
consists of replacing the diagonal, in Nambu space, com- 
ponents of the Green's function with their averages over 
the unit cell of the vortex lattice, while keeping the exact 
spatial dependence of the off-diagonal, Gor'kov compo- 
nents. It was developed for conventional superconductors 
by Brandt, Pesch, and Teword t 84 ' 85 , who showed that 
the replacement is valid since the Fourier components of 
the Green's function (in reciprocal lattice vectors of the 
vortex lattice, K), vary as Gk oc exp(— A 2 K 2 ), where 
A 2 = hc/eB is the magnetic length, which is of order 
of the intervortex distance. The method gives excellent 
agreement with experimental results on both thermody- 
namic and transport properties superconductors near the 
upper critical fiel d 86 ' , and remains semi-quantitatively 
correct in s-wave superconductors down to fields of less 
then half of H c ^^. It fails at the lowest fields, when the 
unpaired electrons are localized in the vortex cores, and 
consequently cannot be described by the propagators av- 
eraged over the (much greater) unit cell size. In this low 
field limit, the method gives and artificially enhanced be- 
havior of the thermal conductivity as it treats the local- 
ized states as extended, and generically produces power 
law increase in k(H, T), while both the expected and the 
experimentally observed initial increase in n(H,T) is ex- 
ponentially small. 

In contrast, the extension of the BPT method to the 
nodal superconductor s 89 ' 90 ' 91 ' 92 ' 93 remains valid down to 
lowest fields since even in that regime the transport is 
dominated by the extended states. It gives results qual- 
itatively and quantitatively consistent with the Doppler 
shift metho d 94 ' 95 , and hence describes the properties over 
nearly the entire range of the temperatures and fields. 
We employ this method here. 

We extend the approach or Refl93l to the two-band 
model. The matrix Green's functions of the electron and 
hole bands are coupled by the self-consistency equation 
on the order parameter and the T-matrix as shown in 
Eq.© and Eqs.([A"8 ]) - (|Allj) . We model the vortex lattice 
as 9 ^ 

MR^)=Y, C ^^y kvV Fo (^-p^j , (16) 

ky 

where F is the ground state oscillator wave function, 
the coefficients Ck y determine the structure of the vortex 
state and the amplitude of the order parameter, and i is 
the band index. Since the bands are treated as separate, 
the thermal conductivity is the sum of the contributions 
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FIG. 9: (Color online) The field dependence of the thermal 
conductivity at T = 0.02T c o for the order parameter with 
nodes, r = 1.3. Left panel: Intraband scattering only for 
different scattering rates near the unitarity limit. Right panel 
top: the same in the presence of strong interband scattering. 
Right panel bottom: evolution of the field dependence with 
increased interband scattering. 



due to the hole and the electron sheets of the Fermi sur- 
face. 

Our approach is most reliable for the states with nodes 
in the gap function, when the results can be trusted over 
essentially the entire field range. Note that even though 
the hole band is always fully gapped in our analysis, at 
low fields the dominant contribution is from the elec- 
tron band with gap nodes. The main feature in the field 
dependence, as shown in Fig. [S] is the pronounced inflec- 
tion point at low fields where a crossover from a rapid 
rise to a slower increase occurs. This result bears strik- 
ing resemblance to the recent measurements on LaFePO 
superconductor—, which were interpreted precisely in the 
framework of the two-band picture, with one band pos- 
sessing nodes in the gap. 

This rapid increase is related to a significant variation 
in the density of states with energy in zero field, shown in 
Fig. [3] The corresponding thermal conductivity for the 
parameter values chosen here also exhibits a shoulder as 
a function of temperature at low T . This shoulder is not 
found experimentally in Ref. l3ll , but we have to keep in 
mind that, if the sample contains non-superconducting 
regions, the residual linear term may not be related to 
the superconducting phase, while the field enhancement 
at low T is still determined by the increase in the number 
of unpaired electrons. 

Note also that when the residual linear- T term is near 
its maximal value in Fig. U which for parameters here 
occurs near r/27rT c o ~ 0.05, the field and temperature 
dependence of k(T, H) is very weak, more reminiscent of 
that of a fully gapped superconductor. The nodes how- 
ever are not lifted up to a higher impurity concentration, 
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FIG. 10: (Color online) The field dependence of the thermal 
conductivity at T = 0.02T c o for the order parameter with 
deep minima, r — 0.9. Left panel: Intraband scattering only. 
Right panel: the strong interband scattering. 



T/2ttT c x) ~ 0.08. The reason for this is clear from com- 
parisons with ReflU: as the slope of the gap near the 
node becomes small, the shape of the gap function devi- 
ates from a simple cosine, and becomes very steep beyond 
the near-nodal region, so that moderate field essentially 
does not excite additional quasiparticlcs. 

In general, we do need to keep in mind that, as the 
nodes are lifted, the applicability of the BPT approxi- 
mation at low fields becomes questionable, but the min- 
imal gap in the regime we show here, T/2ttT c o < 0.15, is 
small, and therefore the method remains reliable to very 
low fields. Quite generally, the energy scales associated 
with the effect of magnetic field on the extended quasi- 
particles are of order E ~ Ay/ H/H c2 , so that for the 
A m i n /A ~ 0.1 we expect the extended states to dom- 
inate the response for H/H C 2 > 0.01. Consequently, 
we trust the approximation even in the regime when the 
nodes arc lifted. 

This argument allows us to extend the treatment to 
the state with the deep minima, rather than the true 
nodes. We consider once again r = 0.9 and show the 
results in Fig. [TO] As can be expected from a probe that 
is sensitive to the amplitude, rather than the phase, of 
the gap function, the overall features are quite similar to 
those for a true nodal gap. The low- field inflection and 
the rapid rise are not as clearly pronounced, consistent 
with the absence of true nodal quasiparticlcs until the 
field is sufficiently high. As pointed out above, for the 
nodeless state the inclusion of strong interband scattering 
leads to a rapid enhancement in the residual density of 
states, and the concomitant increase in the residual linear 
term in the thermal conductivity. 

Finally, in Fig. [11] we show the results for the isotropic 
s± state. While in this case we do not trust our ap- 
proximation at low fields, it is clear that the increase 
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FIG. 11: (Color online) The field dependence of the thermal 
conductivity for the s± state. Left panel: clean case near the 
unitarity limit. Right panel top: influence of the phase shift 
of scattering on the low-field behavior. Right panel bottom: 
evolution of the field dependence with interband scattering. 



in the thermal conductivity with the magnetic field is 
much slower than for the two cases considered above. 
To require an even moderately rapid growth of n(H) at 
low fields requires a substantial residual linear term as 
well as unphysically high interband scattering, see bot- 
tom right panel of Fig. QT] This result strongly suggests 
that the isotropic s± state is incompatible with the re- 
sults of Ref.— on the 122 series of materials, just as the 
results of Refi^i exclude this order parameter structure 
in the LaFePO system. 



VI. CONCLUSIONS 

We have argued that thermal conductivity is the ideal 
probe to resolve current apparent discrepancies between 
various thermodynamic and spectroscopic measurements 
on the Fe-based superconductors. In particular, since 
several of these experiments indicate the existence of low- 
lying quasiparticle states in certain materials, it is impor- 
tant to settle whether or not these excitations extend all 
the way down to the Fermi energy, or whether there is a 
true spectral gap in the system. Thermal conductivity, 
a bulk probe, is currently measurable to lower tempera- 
tures than other probes, so it may be able to settle this 
dispute and also distinguish between two popular sce- 
narios. The two most likely order parameters for these 
systems appear at present to be the isotropic s± state 
proposed by Mazin et al4^, and highly anisotropic A\ g 
states, with nodes or deep gap minima, found in spin 
fluctuation calculations. The former state can be con- 
sistent with the reports of low-lying excitations only if 
pairbrcaking disorder induces an impurity band, while 
the latter are difficult to reconcile with ARPES experi- 
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mcnts indicating a large spectral gap in some materials. 

We have therefore calculated the thermal conductivity 
of a superconductor with A\ g symmetry order parame- 
ter in a model 2-band system, and considered the effects 
of intra- and intcrband disorder. In zero field, we have 
shown that the linear T term which dominates k(T) at 
low temperatures has a coefficient which is nonuniversal 
(unlike the well-known d-wave case) and depends non- 
monotoically on disorder. Details depend on the pre- 
cise order parameter structure of the model pure system, 
and on location of the impurity band in the DOS. The 
linear-T term in zero field found in the LaFcPO mate- 
rial^, is consistent in principle with the linear-T pen- 
etration depth observed in the same system 8 - and sug- 
gestive of nodes in the superconducting order parameter. 
This term is rather large as a fraction of the normal state 
thermal conductivity. While within the present theoret- 
ical approach realistic evaluation of the normal state k 
is difficult, as we have neglected both phonons, which 
would enhance this contribution, and inelastic electronic 
scattering, which would suppress it, it seems likely that 
it would be difficult to account for its size within the cur- 
rent framework, and it is possible, as Ref. [3l| mentions, 
that it is of extrinsic origin. In the 122 systems, the ex- 
tremely small or zero linear-T ter m 28 ' 29 suggests a true 
spectral gap, consistent either with an isotropic s± state 
or a gap with deep minima. 

An examination of the field dependence of the low tem- 
perature thermal conductivity within the BPT approach 
has enabled us to draw further conclusions. The size of 
the initial field dependence seen in experiment rules out 
a clean s± state as a possible candidate for the 122 mate- 
rials. However, as pointed out in the context of other ex- 
periments, pairbrcaking scattering can induce a low-lying 
impurity band in such a state, and produce responses 
similar to highly anisotropic states. We have analyzed 
this situation and found that the amount of pairbreak- 
ing (intraband) scattering required to reproduce the ob- 
served field dependence is large. In most situations this is 
unphysical, both in the sense that the ratio of interband 
to intraband scattering must be tuned to a special value 
(which seems unreasonable in the context of screened 
Coulomb scattering), and because a very large concomi- 
tant T c suppression would be produced. As mentioned 
above, an argument for a sizeable interband scattering 
component may be made for Co-doped systems, but the 
general argument against fine-tuning to a special value 
still holds, and the experimental agreement between the 
behavior of the thermal conductivity on systems with dif- 
ferent dopants indicates the generic features of the ma- 
terial. The evolution of the thermal conductivity with 
Co doping^, however, may be at least in part due to 
the strong interband scattering component. We there- 
fore conclude that the most likely candidate for the or- 
der parameter in the 122 materials is a highly anisotropic 
Ai g state with deep gap minima, probably on the elec- 
tron ("/?") sheets. How this conclusion can be reconciled 
with ARPES experiments is not clear at this writing. We 



emphasize that controlled disorder not associated with 
doping, such as electronic irradiation, would provide the 
best test of the predictions of our theory. 
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APPENDIX A: BASIC FORMALISM 

For low impurity concentrations, one can ignore the 
processes which involve scattering from multiple impurity 
sites. Within this single site approximation, we sum all 
possible scattering events from a single site to calculate 
the disorder-averaged i-matrix, which is then related to 
the one-electron self energy as, 



S(k,w) 



(Al) 



where 



is the impurity concentration, and, for 



isotropic scatterers, Tk,k(w) = T(w), and hence E(u>), 
has no momentum dependence. For a multiband super- 
conductor, scattering from impurities can occur within a 
given band with potential Uu , or between the bands with 
potential Uij, i ^ j. Fig. U2l shows the impurity averaged 
diagrams for the self energy E which occur up to third 
order for a two-band system. Any process which involves 
odd number of interband scatterings does not contribute 
to the self energy, because the Green's function and self- 
energy in the translationally invariant disorder-averaged 
system are diagonal in band index. 

The sum of all the diagrams involving a single impurity 
site can be expressed compactly as 



Ti = 



1 



eff 



(A2) 



where the effective impurity potential for i th band is, 

Uf f = Uu + U i j{Gj) k U ji (A3) 
In a Nambu basis, 



(Gih 



9i,o T o + 9i,i n 



(A4) 
(A5) 



Here g^ a are the Nambu components of the momentum- 
integrated Green's function. The first subscript i(= 1,2) 
in g.j a stands for the band and the second subscript a{= 
0, 1, 2, 3) represents the Nambu channel. 
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Here again the first subscript in £i jQ , represents the band 
index and the second subscript a denotes the Nambu 
channel. 
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APPENDIX B: SPECIAL CASES 
1. Born limit of two band case 



+ 0(U 4 ) 



In the Born limit, we will keep the terms up to second 
order in "£/</' , so the denominator becomes 1 and we get, 



FIG. 12: These are the impurity averaged diagrams, which 
contribute to the self energy of the first band Green's func- 
tion. Here the interband contribution comes through pro- 
cesses, which involve even number of interband scatterings. 
The diagrams also takes into account the order of inter and 
intraband scatterings. Uij is the impurity potential strength, 
wherei, j are the band indexes, i = j denotes the intraband 
and i j denotes the interband potential strength. Gi is the 
bare Green's function. 



The effective potential for the first band may be writ- 
ten 



u ll 



eff 



UnT 3 + U 12 t 3 



Go 



1 



I-U22G2 



T3 



(A6) 



In the above equations, Un,U22 are the intraband im- 
purity potential strengths in band 1 and 2 respectively, 
while U12 is the interband impurity potential strength. 
The i-matrix for, e.g. the first band may now be written 
as 



U e 1 ( S G l 



1 



(A7) 



After some Pauli matrix algebra, we find that the self 
energy components in band and Nambu channels may be 
written as 



£1,0 — 

£1,1 = 

£2,0 = 
V 



^[U? l9lfi + U? 2 g2,o 
0i, 0(^12 - Un U22) 2 (02,0 

[^ii0i,i + E/iatfM 



■ti 



ll)] 



, {U? 2 -Uu U22? (glo-gli)] 



(A8) 
(A9) 



9iA U ': 

£l,2->2,l,a 



= 1-2 Uf 2 (gifi 32,0 - .91,1 ,92,i) (A 

+ {uf 2 - u 11 u 22 ) 2 (.9i,o - ali) (alo - ah) 
- UI2 (ff 2 2 ,o-.9 2 2 i)-^ii (al,o -ah) 



(A10) 
(All) 



£1,0 


= n imp (U^ gi.o + 


U? 2 .92,0) 


£1,1 


= —n-imp (U n g\s 


+ U?2 92,2 


£2,0 


= "Tlimp {U22 52,0 + 


U12 .9l,o) 


£2,1 


= -nimp (U 22 32,2 


+ ^12 91,2 



(Bl) 
(B2) 
(B3) 
(B4) 



2. Strong potential limit 

Due to the presence of an additional band, a new pa- 
rameter comes into play in the unitary limit, the ratio 
of the interband scattering to the intraband scattering. 
This leads to three distinct cases in the large potential 
limit. 

a. Strong potential limit I: Un = U22 > U12 In this 
case, intraband scattering dominates, and in the limit the 
self energies reduces to 



£1,0 

£1,1 
£2,0 
£2,1 



9i.o 



(0i,o - 2 1) 
9i,i 



Hmp , 2 
l.9l,0 



ali) 

02,0 



(9: 



2,0 02, l) 
02,1 



Hmp 



(alo 



ah) 



(B5) 
(B6) 
(B7) 
(B8) 



It is clear that in this limit, states within a given band 
are broadened only by interband scattering processes. 

b. Strong potential limit II: U\\ = U22 = V\ 2 . In 
this very special case, for strong potentials the self ener- 
gies become 



£1/2,0 = ~ n 
£1/2,1 : 



01,0 + 02,0 

P (.91,0 + 02,o) 2 - (01,1 + 02,l) 2 
01,1 +02,1 



(B9) 



imp 



(01,0 +02,o) 2 - (01,1 +02,l) 5 



(BIO) 

(BH) 



In this case, both the bands have identical self energies. 
As discussed in the text, it corresponds to the presence 
of bound states at low energies in the s± state, so we 
devote some attention to it. 



12 



C.C. Tsuei and J.R. Kirtley, Rev. Mod. Phys. 72, 969 
(2000). 

Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 
Am. Chem. Soc. 130, 3296 (2008). 

K. Hashimoto, T. Shibauchi, T. Kato, K. Ikada, R. 
Okazaki, H. Shishido, M. Ishikado, H. Kito, A. Iyo, H. 
Eisaki, S. Shamoto, and Y. Matsuda,Phys. Rev. Lett. 102, 
017002 (2009). 

C. Martin, M. E. Tillman, H. Kim, M. A. Tanatar, S. K. 

Kim, A. Kreyssig, R. T. Gordon, M. D. Vannette, S. Nandi, 

V. G. Kogan, S. L. Bud'ko, P. C. Canfield, A. I. Goldman, 

R. Prozorov, Phys. Rev. Lett. 102, 247002 (2009) 

K. Hashimoto, T. Shibauchi, S. Kasahara, K. Ikada, S. 

Tonegawa, T. Kato, R. Okazaki, C. J. van der Beek, M. 

Konczykowski, H. Takeya, K. Hirata, T. Terashima, Y. 

Matsuda, Phys. Rev. Lett. 102, 207001 (2009). 

R. T. Gordon, N. Ni, C. Martin, M. A. Tanatar, M. D. 

Vannette, H. Kim, G. Samolyuk, J. Schmalian, S. Nandi, 

A. Kreyssig, A. I. Goldman, J. Q. Yan, S. L. Bud'ko, P. 

C. Canfield, R. Prozorov, Phys. Rev. Lett. 102, 127004 
(2009). 

R. T. Gordon, C. Martin, H. Kim, N. Ni, M. A. Tanatar, 
J. Schmalian, I. I. Mazin, S. L. Bud'ko, P. C. Canfield, R. 
Prozorov, Phys. Rev. B 79, 100506 (2009). 
J.D. Fletcher, A. Serafin, L. Malone, J. Analytis, J-H Chu, 
A.S. Erickson, I.R. Fisher, A. Carrington, Phys. Rev. Lett. 
102, 147001 (2009) 

Lin Zhao, Haiyun Liu, Wentao Zhang, Jianqiao Meng, Xi- 
aowen Jia, Guodong Liu, Xiaoli Dong, G. F. Chen, J. L. 
Luo, N. L. Wang, Guiling Wang, Yong Zhou, Yong Zhu, 
Xiaoyang Wang, Zhongxian Zhao, Zuyan Xu, Chuangtian 
Chen, X. J. Zhou, Chin. Phys. Lett. 25, 4402 (2008). 
H. Ding, P. Richard, K. Nakayama, T. Sugawara, 
T. Arakane, Y. Sekiba, A. Takayama, S. Souma, T. Sato, 
T. Takahashi, Z. Wang, X. Dai, Z. Fang, G.F. Chen, 
J.L. Luo, N.L. Wang, Europhys. Lett. 83, 47001 (2008). 
T. Kondo, A.F. Santander-Syro, O. Copie, C. Liu, 
M.E. Tillman, E.D. Mun, J. Schmalian, S.L. Bud'ko, 
M.A. Tanatar, P.C. Canfield, A. Kaminski, Phys. Rev. 
Lett. 101, 147003 (2008). 

D. V. Evtushinsky, D.S. Inosov, V.B. Zabolotnyy, 
A. Koitzsch, M. Knupfer, B. Buchner, G.L. Sun, 
V. Hinkov, A.V. Boris, C.T. Lin, B. Keimer, 
A. Varykhalov, A. A. Kordyuk, S.V. Borisenko, Phys. 
Rev. B 79, 054517 (2009). 

K. Nakayama, T. Sato, P. Richard, Y.-M. Xu, Y. Sekiba, 
S. Souma, G. F. Chen, J. L. Luo, N. L. Wang, H. Ding, T. 
Takahashi, Euro. Phys. Lett. 85, 67002 (2009). 
L. Wray, D. Qian, D. Hsieh, Y. Xia, L. Li, J.G. Checkelsky, 
A. Pasupathy, K.K. Gomes, C.V. Parker, A.V. Fedorov, 

G. F. Chen, J.L. Luo, A. Yazdani, N.P Ong, N.L. Wang, 
M.Z. Hasan, Phys. Rev. B 78, 184508 (2008) 

R. Klingeler, N. Leps, I. Hellmann, A. Popa, C. Hess, 

A. Kondrat, J. Hamann-Borrero, G. Behr, V. Kataev, and 

B. Buechner. larXiv:0808.0708l 

K. Matano, Z. A. Ren, X. L. Dong, L. L. Sun, Z. X. 
Zhao and Guo-qing Zheng, Europhysics Letters 83, 57001 
(2008). 

H. -J. Grafe, D. Paar, G. Lang, N.J. Curro, G. Behr, 
J. Werner, J. Hamann-Borrero, C. Hess, N. Leps, R. Klin- 
geler, and B. Buchner, Phys. Rev. Lett. 101, 047003 



(2008). 

18 K. Ahilan, F.L. Ning, T. Imai, A.S. Sefat, R. Jin, 
M.A. McGuire, B.C. Sales, D. Mandrus, Phys. Rev. B 78, 
100501 (R) (2008). 

19 T.Y. Nakai et al., J. Phys. Soc. Jpn. 77, 073701 (2008). 

20 M. Yashima, H. Nishimura, H. Mukuda, Y. Kitaoka, K. 
Miyazawa, P. M. Shirage, K. Kiho, H. Kito, H. Eisaki, and 
A. Ivo, larXiv:0905.1896l f unpublished) . 

21 L. Shan, Y. Wang, X. Zhu, G. Mu, L. Fang, C. Ren and 
H.-H. Wen, Europhys. Lett. 83, 57004 (2008). 

22 T.Y. Chien, Z. Tesanovic, R.H. Liu, X.H. Chen, and C.L. 
Chien, Nature 453, 1224 (2008). 

23 D. Daghero, M. Tortello, R.S. Gonnelli, V.A. Stepanov, 
N.D. Zhigadlo and J. Karpinski, Journal of superconduc- 
tivity and novel magnetism, 22,553 (2009). 

24 R.S. Gonnelli, D. Daghero, M. Tortello, G.A. Ummarino, 
V.A. Stepanov, J.S. Kim and R.K. Kremer, Phys. Rev. B 
79, 184526 (2009) 

25 Y. Matsuda, K. Izawa, and I. Vekhter, J. Phys.: Cond. 
Mat. 18, R705-R752 (2006). 

26 H. Shakeripour, C. Petrovic, Louis Taillefer, New Journal 
of Physics 11, 055065 (2009). 

27 J. G. Checkelsk y, Lu Li, G. F. Ch en, J. L. Luo, N. L. Wang 
and N. P. Ong, larXiv:0811.4668l 

28 X. G. Luo, M. A. Tanatar, J.-Ph. Reid, H. Shakeripour, 
N. Doiron-Leyraud, N. Ni, S. L. Budko, P. C. Canfield, 
Huiqian Luo, Zhaosheng Wang, Hai-Hu Wen, Ruslan Pro- 
zorov, and Louis Taillefer, arXiv:09 04.4049l 

29 L. Ding, J. K. Dong, S. Y. Zhou, T. Y. Guan, X. Qiu, C. 
Zhang, L. J. Li, X . Lin, G. H. Cao, Z. A. Xu, and S. Y. Li, 
larXiv:0906.0138l 

30 Y. Machida, K. Tomokuni, T. Isono, K. Izawa, Y. Naka- 
jima, and T. Tamegai, larXiv:0906.0508l 

31 M. Yamashita, N. Nakata, Y. Senshu, S. Tonegawa, K. 
Ikada, K. Hashimoto, H. Sugawara, T. Shibauchi, and Y. 
Matsuda, larXiv:0906.0622l 

32 Note that there is substantial disagreement in the 
low-temperature behavior of Co-doped samples between 
Refs.[l]and[37|. 

33 G. Preosti and P. Muzikar, Phys. Rev. B 54, 3489 (1996). 

34 A. A. Golubov and I.I. Mazin, Phys. Rev. B 55, 15146 
(1997). 

35 Y. Senga and H. Kontani, J. Phys. Soc. Jpn. 77, 113710 

(2008) . ;New J. Phys. 11, 035005(2009). 

36 Y. Bang, H.-Y. Choi, and H. Won, Phys. Rev. B 79, 054529 

(2009) . 

37 M. A. Tanatar, J. P. Reid, H. Shakeripour, X. G. Luo, N. 
Doiron-Leyraud, N. Ni, S. L. Bud'ko, P. C. Canfield, R. 
Prozorov, Louis Taillefer, arXiv:0907.1276 (unpublished) 

38 V. Mishra, G. Boyd, S. Graser, T. Maier, P.J. Hirschfeld, 
and D.J. Scalapino, Phys. Rev. B 79, 094512 (2009) 

39 K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. Kon- 
tani and H. Aoki, Phys. Rev. Lett. 101, 087004 (2008). 

40 F. Wang, H. Zhai, Y. Ran, A. Vishwanath and D.-H. Lee, 
Phys. Rev. Lett. 102, 047005 (2009). 

41 S. Graser, T. A. Maier, P. J. Hirschfeld, and D. J. 
Scalapino, New J. Phys. 11, 025016 (2009). 

42 T. Maier, S. Graser, P.J. Hirschfeld, and D.J. Scalapino, 
Phys. Rev. B 79, 224510 (2009). 

43 A.V. Chubukov, M.G. Vavilov, A.B. Vorontsov, 
larXiv:0903.5547l 



13 



44 R. Thomale, C. Piatt, J. Hu, C. Honerkamp, B. A.Bernevig 
larXiv:0906.4475l 

45 I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008). 

46 A.V. Chubukov, D. Efremov, and I. Eremin, Phys. Rev. B 
78, 134512 (2008). 

47 V. Ambegaokar and A. Griffin, Phys. Rev. 137, A1151 - 
A1167 (1965). 

48 C. Cao, P.J. Hirschfeld, H.-P. Cheng, Phys. Rev. B 77, 
220506(R) (2008). 

49 PA. Lee, Phys. Rev. Lett. 71, 1887 (1993). 

50 M. Graf J.A. Sauls, Phys. Rev. Lett. 53, 15147 (1996). 

51 M. Norman and P.J. Hirschfeld, Phys. Rev. B 53, 5706 
(1996). 

52 M. T. Beal-Monod and K. Maki, Physica C 265, 309 

(1996) . 

53 L. S. Borkowski and P.J. Hirschfeld, Phys. Rev. B 49, 
15404 (1994). 

54 R. Fehrenbacher and M. R. Norman, Phys. Rev. B 50, 
3495 (1994). 

55 A.B. Vorontsov, M.G. Vavilov, and A.V. Chubukov, Phys. 
Rev. B 79, 140507 (2009). 

56 D. Parker, O.V. Dolgov, M.M. Korshunov, A. A. Gohibov, 
and LI. Mazin, Phys. Rev. B 78, 134524 (2008). 

57 A.V. Chubukov, D. Efremov, and I. Eremin, Phys. Rev. B 
78, 134512 (2008). 

58 D. Markowitz and L. P. Kadanoff, Phys. Rev. 131, 563 
(1963). 

59 A. A. Abrikosov and L.P. Gor'kov, Zh. Eksp. Teor. Fiz 39, 
1781(1960) [Sov. Phys. JETP 12, 1243(1961)]. 

60 L.P. Gorkov and PA. Kalugin, Pisma Zh. Eksp. Teor. Fiz. 
41, 208 (1985) [JETP Lett. 41, 263 (1985)]. 

61 K. Ueda and M. Rice, 1985, in Theory of Heavy Fermions 
and Valence Fluctuations, edited by T. Kasuya and T. Saso 
(Springer- Verlag, Berlin), p. 216. 

62 M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 

63 S. Skalski, O. Betbeder-Matibet, and P. R. Weiss, Phys. 
Rev. 136, A1500 (1964). 

64 P.J. Hirschfeld, P. Woffle and D. Einzel, Phys. Rev. B 37, 
83 (1988). 

65 We have averaged the xx and yy components to simulate 
the ferropnictide with order parameters on /3i and P2 sheets 
rotated by 7r/2 in local coordinates with respect to one 
another. 

66 The approximate factor of 2 difference between the kink 
temperature and the antinodal order parameter energy en- 
ergy arises from the width of the thermal quasiparticle dis- 
tribution. 

67 F. Gross, B.S. Chandrasekhar, D. Einzel, P.J. Hirschfeld, 
K. Andres, H. R. Ott, J. Beuers, Z. Fisk and J.L. Smith, 
Z. Physik B 64, 175 (1986). 

68 I. Koztin and A.J. Leggett, Phys. Rev. Lett. 79, 135 

(1997) . 

69 Y. Kamihara, et al. J. Amer. Chem. Soc. 128, 10012 



(2006). 

70 A.I. Coldea et al., Phys. Rev. Lett. 101, 216402 (2008). 

71 P.J. Hirschfeld and N. Goldenfeld, Phys. Rev. B. 48, 
4219(1993). 

72 L. S. Borkowski, P.J. Hirschfeld, and W.O. Putikka, Phys. 
Rev. B 52, 3856 (1995). 

73 R. Movshovich, M. Jaime, J. D. Thompson, C. Petrovic, 
Z. Fisk, P. G. Pagliuso, and J. L. Sarrao , Phys. Rev. Lett. 
86, 5152 (2001) 

74 M. Sutherland, D. G. Hawthorn, R. W. Hill, F. Ronning, 
S. Wakimoto, H. Zhang, C. Proust, Etienne Boaknin, C. 
Lupien, Louis Taillefer, Ruixing Liang, D. A. Bonn, W. 
N. Hardy, Robert Gagnon, N. E. Hussey, T. Kimura, M. 
Nohara, and H. Takagi, Phys. Rev. B 67, 174520 (2003). 

75 E. Boaknin, M. A. Tanatar, Johnpierre Paglione, D. 
Hawthorn, F. Ronning, R. W. Hill, M. Sutherland, L. 
Taillefer, J. Sonier, S. M. Hayden, and J. W. Brill, Phys. 
Rev. Lett. 90, 117003 (2003) 

76 Louis Taillefer, Benoit Lussier, Robert Gagnon, Kamran 
Behnia, and Herv Aubin, Phys. Rev. Lett. 79, 483 (1997) 

77 G. E. Volovik, JETP Letters 58, 469 (1993). 

78 C. Kiibert and P. J. Hirschfeld, Solid State Commun. 105, 
459 (1998). 

79 I. Vekhter, P. J. Hirschfeld, and E. J. Nicol, Phys. Rev. B 
64, 064513 (2001). 

80 C. Kiibert and P. J. Hirschfeld, Phys. Rev. Lett. 80, 4963 

(1998) . 

81 M. Franz, Phys. Rev. Lett. 82, 1760 (1999). 

82 I. Vekhter and P. J. Hirschfeld, Physica C 341-348, 1947 
(2000). 

83 I. Vekhter, L. N. Bulaevskii, A. E. Koshelev, and M. P. 
Maley, Phys. Rev. Lett. 84, 1296 (2000). 

84 U. Brandt, W. Pesch, and L. Tewordt, Z. Phys. 201, 209 
(1967). 

85 W. Pesch, Z. Phys. B 21, 263 (1975). 

86 A. Houghton and K. Maki, Phys. Rev. B 4, 843 (1971). 

87 P. Klimesch and W. Pesch, J. Low Temp. Phys. 32, 869 
(1978). 

88 E. H. Brandt, J. Low Temp. Phys. 24, 409 (1976). 

89 I. Vekhter and A. Houghton, Phys. Rev. Lett. 83, 4626 

(1999) . 

90 H. Kusunose, Phys. Rev. B 70, 054509 (2004). 

91 A. B. Vorontsov and I. Vekhter, Phys. Rev. Lett. 96, 
237001 (2006). 

92 A. B. Vorontsov and I. Vekhter, Phys. Rev. B 75, 224501 
(2007) 

93 A. B. Vorontsov and I. Vekhter, Phys. Rev. B 75, 224502 
(2007) 

94 T. Dahm, S. Graser, C. Iniotakis, and N. Schopohl, Phys. 
Rev. B 66, 144515 (2002). 

95 G. R. Boyd, P. J. Hirschfeld, I. Vekhter, and A. B. 
Vorontsov, Phys. Rev. B 79, 064525 (2009) 



